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We consider a singly quantized off-axis straight vortex in a rotating dipolar ultracold gas in 

tiie Thomas-Fermi (TF) regime. We derive analytic results for small displacements and perform 

numerical calculations for large displacement within the TF regime. We show that the dipolar 

interaction energy increases (decreases) as the vortex moves from the trap center to the edge in an 

("^ ■ oblate (a prolate) trap. We find that for an oblate (a prolate) trap, the effect of the dipole-dipole 

fSJ ' interaction is to lower (raise) both the precession velocity of an off-center straight vortex line and 

K.^ , the angular velocity representing the onset of met ast ability. 
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I. INTRODUCTION 



__ 'I The first experimental detection of a vortex in a dilute alkali- atomic gas Bose-Einstein condensate (BEC) was made 

d • by Matthews et al. in 1999 using ®^i?6 atoms [1] and theoretical predictions on the main features of the vortex states 
(^' have been shown to agree with experiments 2] (references therein) . Rb has a small dipole moment while chromium 
(.J , atoms posses a larger permanent magnetic dipole moment, which leads to significant dipolar interactions in addition to 
^ the usual short-range interactions. The successful Bose-Einstein condensation of ^^Cr atoms has stimulated a growing 

^ , interest in the study of BEC with nonlocal dipole-dipole interactions [3-5]. This nonlocal character has remarkable 
consequences for the physics of rotating dipolar gases |6|-|20|. It has been shown that the critical angular frequency 
for vortex creation may be significantly affected by the dipolar interaction [y|. In addition, dipolar gases under fast 
"j^ . rotation develop vortex lattices, which due to the dipolar interaction may be severely distorted 11 2| . and even may 
Ch ' change its configuration from the usual triangular Abrikosov lattice into other arrangements |13l . Il4| . It was shown 
I , that, the dipolar interaction may significantly modify the vortex line stability. Under appropriate conditions, the 
' O ■ dispersion law for transverse modes shows a rotonlike minimum, which for sufficiently large dipolar interaction may 
G ' reach zero energy, destabilizing the Kelvin waves 15:]. In the TF limit, dipole-dipole interaction changes the stability 
and the instability conditions and the possibility of vortex lattice formation for a rotating dipolar BEC in an elliptical 
trap HI]. 

The long-range and anisotropic interactions introduce rich physical effects, as well as opportunities to control BECs. 
In a prolate dipolar gas with the dipoles polarized along the z-axis, the dipolar interaction is attractive, whereas it 
is repulsive for an oblate dipolar gas. As a result, the sign of the dipolar mean-field energy can be controlled via the 
trap aspect ratio. In this paper, we will consider an off-axis vortex lines in an oblate dipolar BEC with the dipoles 
r^ , aligned in the z direction by an external field. In an oblate condensate, vortex line can be approximated as straight. 
fT^ ■ This is not the case for a prolate condensate. In that case, vortex lines are twisted. In the case of short range contact 
, -/ ] interaction, an off- axis vortex in a BEC was studied in detail |2l| - t26| . We will analyze the effects of the dipolar 
(«— ^ ■ interaction on the physics of an off-axis vortex. Specifically, we shall assume that a straight singly quantized vortex 
f~^ line in an oblate trap is displaced from the trap center of the Cr condensate with transverse coordinates xq and j/o- 
T-H This paper is structured as follows. Sec. II reviews the TF solution for a dipolar condensate. Section III investigates 

s^ ' the singly quantized straight vortex line in the presence of the dipole-dipole interaction. The last section discusses 
the results. 



II. TF SOLUTION 

In this section, we will review the TF solution for a dipolar gas. Interparticle interaction potential in dipolar 
gases includes both a short-range Van der Waals and a long-range dipole-dipole terms. Because of the long-range 
character of the dipole-dipole interaction, scattering properties at low energies are significantly changed. In the case 
of a short-range interaction, only the s-wave scattering is important at low energies. However, in the case of a long- 
range interaction, all partial waves contribute to scattering. Within the mean-field description of the condensate, the 
interaction potential is well described by the following model potential J27l - l31| 

V = gS{r) + ^il~3cos^e) , (1) 



where o = , a^ is the scattering length, d is the electric dipole moment (the results are equally valid for magnetic 

m 
dipoles), r is the vector connecting two dipolar particles and 6 is the angle between r and the dipole orientation. In 

this study, we will suppose that the dipoles are polarized along the z-axis. 

For dipolar condensates, it is useful to introduce a dimensionless parameter that measures the relative strength of the 

dipolar and s-wave interactions 

--t (2) 

where the coupling Cdd = MoA*^- Chromium atoms posses an anomalously large magnetic dipole moment ficr = 6/iB 
{hb is the Bohr magneton), while Rb has a dipole moment equals to fiub = 1a*s [31]. It has been shown in :32] that 
in the Thomas- Fermi limit a dipolar BEC also is stable as long as < Edd < 1- 

Consider a dipolar BEC of N particles with mass m and electric dipole d oriented in the z-direction by a sufficiently 
large external field. At sufficiently low temperatures, the description of the ground state of the condensate is provided 
by the solution of the Gross-Pitaevskii (GP) equation 

g\^{r,t)\^ + d^^dd{r)]'i'{r,t) , (3) 



dt ~ 


\ 2m 


where Vr is the trap potential 






m , 

Vrp = — CJ 



(P^+7V), (4) 

where p^ — x^ + y^ and 7 = —^ is the trap aspect ratio and ^ddM — / ^^r' — j — |\I'(r',i)p is the mean-field 

wx J |r-r'|'^ 

potential due to dipole-dipole interactions. 

The equation ([3]) is an integro-differential equation since it has both integrals and derivatives of an unknown wave 

function. This equation can be solved analytically if we assume that the zero-point kinetic energy associated with 

the density variation becomes negligible in comparison to both the trap energy and the interaction energy between 

atoms. In this case, the kinetic energy term can be omitted in the equation. This approximation is known as TF 

approximation. Eberlein et al. showed that a parabolic density profile remains an exact solution for an harmonically 

trapped vortex- free dipolar condensate in the TF limit |32| . 

ribgir) = no{l - ^ - j^) (5) 

In the absence of dipolar interaction, the condensate aspect ratio, k = — , and the trap aspect ratio, 7, match. However, 

1j 

the presence of the dipolar interaction changes the condensate aspect ratio. It satisfies the following equation 

S^^'sddii^ + l)P% - 1] + [Sdd - l){^^ - 7') - (6) 

where f{hi) for oblate case, {k> 1), is given by 



.2(/\ a arctan -J k^- 



2 + k2(4-6 



In the TF regime, the mean-field potential integral due to dipole-dipole interactions, ^^^(r), can be evaluated in the 
spheroidal coordinates [32|. The result can be expressed in cylindrical coordinates. 

Let us now write the energy expression for a dipolar gas. The total energy for a dipolar gas can be written as 

Etot — Ekin + Etrap + Egw + Edd , (9) 

where Ekin is the kinetic energy 

Ek^n = -^j d\^*{v)y^^{v) , (10) 



Etrap is the trap energy 

Etrap= I d\\^{v)\''VT , (11) 

Esw is the energy due to the short range interaction 

Es^^^-jd\\^{v)W (12) 

and Efij^ is the energy due to the long range dipole-dipole interaction 

Em = \ j d?rSr'\^{v)\^UM{v -v')\^{v')\^ = \ j d?r n(r)$dd(r) . (13) 

In the TF approximation, the kinetic energy term ([TU| is neglected. The total energy associated with the vortex-free 
Thomas-Fermi solution is given by [32j] 

E,,, = -mulR\2 + ^) + — -^ (1 - e,,f{^)) (14) 

III. DIPOLAR CONDENSATE WITH A VORTEX 

Consider a single straight vortex line at a position po along the z-axis. In this case, the wave function, normalized 
to the total number of atoms / |\I'|^d'^r = A^, is given by 



*(p, (j,, z) = v/^ e'^^P^P"^ , (15) 

where n(r) = |^(r)p is the density. The expression of the phase S{p,po) characterizing the circulating flow around 
the vortex line is given by |2| 

S{p,po) — arctan I ) (16) 

where Pq — Xq + i/q. The corresponding irrotational flow velocity is given by 

v=l\/S{p,po) (17) 

m 

There is a singularity on the vortex line, p = po, where the velocity diverges. However, the particle current density, 
J = nv, vanishes as p — )► po- When a quantized vortex is present at the position po, the density drops to zero at 
the center of the vortex core whose size is determined by the parameter [3. For a centered vortex in a BEC without 
dipole-dipole interaction, the parameter /? is given by 

where d± = y'Ti/muj^ is the mean oscillator strength. The TF limit holds when R is large compared to d±. The TF 
length scale reads /3 << d± << R. The vortex core size increases with po- The parameter Pswifo) characterizing the 
small vortex core at the position po in a BEC without the dipole dipole interaction is [22| 

Having written the expression of the phase, for a straight off-center vortex, let us now find the density n(r). The 
repulsive interactions and the repulsive dipolar interaction (for oblate case) significantly expand the condensate, so 
that the kinetic energy associated with the density variation becomes negligible compared to the trap energy and 



interaction energies. In the TF regime, the density profile of a condensate with a straight off-axis vortex line at po is 
given by lU 

where n(r) — when the right hand side is negative and /3, R, and L are variational parameters that describe the size 
of the vortex core, and the radial and the axial sizes, respectively. Note that the density function ((20)) behaves like 

\p — poP//^^ when p « (3 and like (1 — — — —) when p » (3. These parameters will be calculated by minimizing 

IX Ij 

the energy functional. The central density no can be found using the normalization condition 

^0 = y^ (eOi?^ + 25^2 (^4i?2(3 ln(|) - 4) + 9pg(3 - 2 ln(|))) ) (21) 

We don't include an image vortex because the form of the TF condensate density ensures that the particle current 

density vanishes at the surface. 

Let the total angular momentum for a singly quantized vortex line along the trap axis at the position po be L^, 

{Lz = m rvij,n(r)d r). Then the corresponding energy of the system in the rotating frame is E' = E — ^L^, where 

E is the energy in the non-rotating frame. If we denote the energy of the BEC in its ground state without a vortex 
by Eq and the extra energy needed to generate a vortex by AE, then E = Eq + AE. We can now write the energy 
of the vortex state in the rotating frame as 

E' = Eo + AE- flL^ (22) 

A vortex is generated if E' is smaller than Eq. In other words, a vortex is formed above a certain critical value of the 
rotation frequency. The critical rotational velocity is given by 

fie = ^ (23) 

It should be noted that a vortex lattice starts to appear when the rotation frequency is further increased. 

We proceed by minimizing the total energy with respect to the three variational parameters, R, k and /3. Let us now 

calculate the kinetic, trap, s-wave and dipole-dipole interactions energies separately. Since (3 is small, we will neglect 

terms of order /3 and higher. The energy integral can be evaluated analytically up to the second order of po- Below 

we will obtain analytical expression for small pg. In the following section, we will perform numerical calculations for 

large po in the TF limit. 

Let us firstly obtain the kinetic energy 

n^nnoR f ^„ , ^„,^ , o52m__.2, ^^-^ ( o«2m„/2x\^2 



Ehu = g^^ ( -22 + 12(1 + 3/3^) In(^) - 27/3^ + 18 ( 2 - (1 + 2/3^) In(^) ) p^ ) (24) 



where we have defined 



^ = |, p-o = f . (25) 

Using the expression (lll[) . the trapping energy is straightforwardly evaluated to be 

E^ra, = ^^^^^5^ (y(2.^ + 7^) - g(15.^ - 7^(23 - 151n(|))) + ^ {2%.^ - U^' + 61n(|)(7^ - 2^^)) pi 



In the similar way, the formula (|12p yields the s-wave interaction energy 



(26) 



E.. . 5^ n + i^^= - V' ,„,|) - ? I 25 - 12 l„,i) I i>i I (27) 



15k V7 15 ^ 6 V P 

Let us now calculate the dipole-dipole interaction energy. Since /3 is small, the dipolar energy function can be 
approximated as [6| 



Em ^\j d\n,g{v)^\%v)+ j Srn^iv)^'^ 



;(r) (28) 



where ribg was defined in ([5]) and riy is defined as 

Note that n(r) = ri6g(r) + n,j(r). Hence, the dipolar interaction energy becomes Edd — E^dd + ^Edd, where 



^^^|g^ (-f /(^) + ^^ (601n(|) - 122 + -g^(62 - 245.^ + 30(5«^ - 2) hr(|) 



l\Edd = ;:^ 50- 24hi(-^) + ^r -(6 + 69k -36k ln(-)) Po (31) 

yK \ p K^ — i P / 

We have obtained the energy expressions up to the second order of fractional vortex displacement, po- Note that the 

central density no in these expressions also includes pg pTj) . Up to the order of /3^, they agree with the results p in 

the limit po^ 0. 

As can be seen, the kinetic energy decreases with po- The kinetic energy goes to zero as pq^R since TF density 

vanishes at the surface. Note that the description of a vortex close to the boundary is outside the scope of the present 

approach, since TF approach doesn't work close to the surface. The dipole-dipole interaction increases with po for an 

oblate trap while decreases with po for a prolate trap. The kinetic energy depends on pg while the dipolar, trap and 

the s-wave interaction energies depend on /3 pg. 

Before embarking on a specific example, let us study the energy expressions qualitatively for an oblate trap. Firstly, 

let us investigate roughly how the total energy is distributed among kinetic, dipolar, trap and s-wave interaction 

energies. The ratio between the kinetic energy and the trap energy is of order /3 ; Ekin^ P Etrap- The trap and 

the s-wave interaction energies are comparable to each other; Esw~ Etrap- The ratio between dipolar and the s-wave 

interaction energies is of order edd] Edd~ ^ddEsw 

Secondly, let us investigate how the excess energy /S.E needed to generate a vortex is distributed. Consider first a 

central vortex, pg = 0. The excess energy for the trap, dipolar and s-wave interaction energies vary as /3 . However the 

excess kinetic energy is of order the kinetic energy, AEkin^Ekin^ ^"^Etrap- Hence, AEkin^ AEtrap- We emphasize 

that the excess energy for dipole-dipole and s-wave interaction energies are negative. Hence, the effect of increasing 

dipole moment and scattering length is to decrease the critical angular velocity flc- The relations between the excess 

energies for dipolar, trap and s-wave terms are given by AEdd^ ^dd AEsw and AEtrap~ — AEsw 

Finally, let us mention how the energy changes with the position of an off-axis vortex, pg. The kinetic energy and the 

trap energy decrease with pg while the dipole-dipole and the s-wave interaction energies increase. Furthermore, the 

total energy decreases with pg. 

Let us briefly study qualitatively the energy expressions for a prolate trap with a straight vortex line. If we neglect 

vortex bending effect, we can use the above energy expressions. In this case, the dipolar interaction energy is negative. 

However, the excess energy for the dipolar interaction is positive. Hence, the effect of increasing dipole moment for a 

prolate trap is to increase the critical angular velocity J7c. Finally, the dipole-dipole interaction energy decreases with 

increasing pg. As a result, the effect of dipolar interaction is to repel an off axis vortex away from the trap center for 

a prolate trap while attract it to the trap center for an oblate trap. 

In what follows, we will give an explicit example for a straight off-axis vortex for an oblate trap. 

IV. RESULTS 

We will study a dipolar BEG with a single vortex in an oblate trap including 150000 ^^Cr atoms. We take the 
numerical values used in the reference [g| to compare the off-center vortex to the central vortex. The trap frequencies 
are lli± = 27r x 200 rad/s and w^ = 27r x 1000 rad/s for 7 = 5. The harmonic oscillator length of the trap along the 
radial direction is d± = 0.986pTO. The magnitude of the magnetic dipole interaction for Cr is Cdd = Mo(6M-b)^- For 
small values of pg, we will use the analytical results obtained in the previous section. For large values pg, numerical 
computation within the TF limit will be performed. 

Let us firstly analyze three variational parameters, (3, k and R. The density of the condensate drops to zero at the 
center of the vortex core whose size is equal to /?. It is very small compared to the radial size of the condensate. The 
smallness of P ensures that the vortex affects the density only the immediate vicinity of the core. Fig-[T] depicts the 
fractional vortex core size, /? = /?/i?, versus the scattering length. The solid curve corresponds to a central vortex, 
while the dashed curve to an off-center vortex with pg = 0.4_R. The parameter /3 is bigger in the presence of an off-axis 
vortex. As can be seen from the figure, the fractional vortex core size decreases with increasing scattering length. 



This can be understood as follows. The radial size increases as scattering length is enlarged. The vortex core size is 
inversely proportional to radial size. So, we conclude that /3 decreases with a^. 

Similarly, Fig-[2] show the aspect ratio k (left) and the radial size (right) of the condensate R versus the scattering 
length, respectively for po = and po — 0.4i?. Contrary to the case of vortex core size /3, the parameters k and R 
don't change appreciably with po when as > 50ao, where ao is Bohr radius. Hence, the curves lie on top of each other 
in Fig-H 

Having discussed the three variational parameters, let us now study the total energy of a dipolar condensate. Figures 
13] and m show the total energy as a function of vortex position for fixed a^ = lOOao and scattering length for fixed 
Po = 0.2 in a non-rotating oblate trap (f2 = 0), respectively. The solid curve corresponds to a condensate with s-wave 
plus dipolar interactions while the dashed line corresponds to a condensate with pure s-wave interaction. The total 
energy is bigger when £^^7^0. This is because the dipole-dipole interaction energy is positive in an oblate trap. The 
total energy of the system attains a maximum when po = for both cases. As the off-axis vortex moves to the edge 
of the condensate, the total energy decreases. More specifically, the kinetic and trap energies decrease with po while 
the dipolar and s-wave interaction energies increase with po. In fact, higher than a specific value of edd, dipolar 
interaction becomes more dominant, so total energy increases with po. We calculate that this happens when edd > 1. 
As mentioned in [6], however, the condensate enters an instability region when edd > 1. As can be seen from the figure 
m the energy differences between the two cases decreases when the scattering length is increased. This is because edd 
is decreased with increasing scattering length ([2]). 

For the investigation of the vortex generation, not the total energy but the excess energy AE associated with the 
presence of an off- axis straight vortex is more important. Fig-[S] compares the excess energy of the condensates with 
Edd = 0.15 (solid curves) and edd = (dashed curves) as a function of a fractional vortex displacement. Different 
curves represent different fixed values of the external angular velocity J7. The top solid and dashed curves correspond 
to ri = 0, where fl increases as one moves towards the lowest curve with fl = 0.08 o;^ and Q, — 0.140 lu^. Note that 
the critical rotation frequency is Vtc = 0.124 uj± (i7c = 0.119 lo^) when edd = [edd = 0.15). As can be seen from 
the figure, the dipolar interaction lowers /S.E compared to the pure contact interaction. It is of great importance 
to note that although dipolar interaction is positive for an oblate trap, the excess dipolar energy is negative. As 
Po is increased, the curves for edd — 0.15 and edd — start to coincide. The top two curves show that the excess 
energy /S.E decreases monotonically with increasing po, with negative curvature at po — 0. So, a central vortex is 
unstable to infinitesimal displacements. The presence of dissipation will move an off-axis vortex toward the edge of 
the condensate. If the trap is rotated with angular velocity, 51, then the energy of a vortex decreases. Inspection of 
Fig. [S] reveals that with increasing rotation speed, the function /S.E flattens. At a special value of rotation frequency, 
rim, curvature of the function AE' becomes zero at po = Q. Hence, above an angular velocity $7^, the vortex attains 
a local minimum. The central position is not globally stable but locally stable. One of the results of this paper is 
that the presence of dipolar interaction lowers Vim for an oblate trap. Let us look at the lowest curves in Fig. [5j In 
this case, appearance of a vortex becomes energetically favorable since /S.E < 0. The central vortex is both locally 
and globally stable relative to the vortex-free state. A vortex initially placed off-center will follow a path of constant 
energy under the action of the Magnus force, which is proportional to the gradient of the energy in the radial direction. 
The precession velocity of a displaced vortex of a nonrotating trap increases with the vortex displacement. Hence, 
a vortex near the surface precesses more rapidly than one near the center. Another result of this paper is that the 
precession velocity of a displaced vortex is lowered in the presence of the dipolar interaction in an oblate trap. On 
the contrary, it is raised in a prolate trap (ignoring vortex bending effect). Note that the precession velocity around 

FiE f)E/fin 
the center for a nonrotating trap, w, can be calculated using w = — — = — — , where E is the energy and L is the 

oL oL/opo 
angular momentum |2|. |22|. For a condensate in rotational equilibrium at angular velocity 51, the original precession 
frequency is altered to a; —> a; — 51 [2] . 

Finally, in Fig-[6l we have examined the critical angular velocity of the condensate for 7 = 5 and 7=10. The critical 
angular velocity above which a vortex state is energetically favorable depends on 7. As can be seen, 51c increases with 
decreasing 7. For stirring frequencies below 51c, no vortex can be nucleated. The presence of dipole dipole interaction 
decreases 51c for an oblate trap. 

We have found that flc, 51^ and precession velocity decrease (increase) in an oblate (a prolate) trap. 
This can be understood simply as follows. The dipolar mean field potential has a parabolic profile 
$1^9 (p, z) = noCdd/3L^(1.21 - 0.04p^ - 1.83z^) when edd = 0.15 and 7 = 5®. This potential has the same inverted 
parabola shape as in the case of contact interactions. So, we conclude that there is a similarity between dipolar and 
non-dipolar BEG in the TF regime. The difference is in the expressions for the radial and axial size. It is well known 
that the contact interaction with positive scattering length decreases the critical angular frequency 51c (^c = w_l for 
a noninteracting trapped gas). In the similar way, n„i/^± and precession velocity decrease with increasing scattering 
length. So, we conclude that inclusion of dipolar interaction in an oblate trap reduces 51c, ^m and precession velocity 
in the TF regime. Furthermore, if we ignore vortex bending effect, the mean field dipolar potential for a prolate 



trap has the same form as the mean field contact potential with negative scattering length. In contrast to the case 
for repulsive interactions, ficj ^m smd precession velocity increase in the presence of attractive contact interactions. 
Analogously, we conclude that dipolar interactions increase them in a prolate trap. 

In this paper, an ofF-axis vortex lines in an oblate ^^Cr BEC polarized along the z direction have been studied. The 
effects of the dipolar interaction on the physics of an ofF-axis vortex have been analyzed. It was shown that the 
condensate aspect ratio, k, and the radial size, i?, remain almost the same in the presence of an off-axis vortex when 
Os > 50ao. The dipolar interaction raises (lowers) the total energy in an oblate (a prolate) trap. On the contrary, the 
excess dipolar energy needed to generate a vortex decreases (increases) in an oblate (a prolate) trap. It was found that 
the angular velocity 51™ representing the onset of metastability and the critical angular velocity of the condensate ilc 
are lowered (raised) in an oblate (a prolate) trap. Finally, it was proven that the effect of the dipole-dipole interaction 
is to lower (raise) the precession velocity of an off-axis straight vortex line around the center in an oblate (a prolate) 
trap. 
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FIG. 1: The fractional vortex core size for a cental (p = 0) and an off-axis (p = 0.4) vortices versus the scattering length for a 
dipolar BEG with a vortex in an oblate trap with 7 = 5. The scattering length is measured in units of Bohr radius, ao- 
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FIG. 2; For p = Q and p = 0.4, the aspect ratio (left) and the radial size (right) of a dipolar BEG with a vortex in an oblate 
trap with 7 = 5 as a function of the scattering length. The scattering length is measured in units of Bohr radius, ao. 
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FIG. 3: The total energy of a non-rotating dipolar BEG in an oblate trap with 7 = 5 for the fixed scattering length a^ = lOOao 
as a function of vortex displacement. The solid curve indicates both dipolar and s-wave interaction with tdd ~ 0.15 while the 
dashed curve indicates only s-wave interaction. 
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FIG. 4: The total energy of a non-rotating dipolar BEG with a vortex in an oblate trap (7 = 5) with po = 0.2R as a function 
of the scattering length (in units of Bohr radius). The solid (dashed) curve is for the condensate with both dipolar and s-wave 
interaction (only s-wave interaction). 




FIG. 5: The increased energy AE in units of Nhuj± in the rotating frame associated with the presence of an off-axis straight 
vortex as a function of a fractional vortex displacement in an oblate trap. The solid (dashed) curves correspond to edd ~ 0.15 
(edd = 0). Different curves represent different fixed values of the external angular velocity Q. The top solid and dashed curves 
corresponds to il = 0, where fl increases as one moves towards the lowest curve with il = 0.08 u)± fl = 0.14 ijj±. 
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FIG. 6: The critical angular velocity of a condensate with a vortex for 7 = 5 and 7 = 10 as a function of the scattering length. 
Qc is measured in units of ijj± 



